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ABSTRACT 

We analyze the stability of a dilute plasma with thermal and composition gradients in the limit where con- 
duction is slow compared to the dynamical timescale. We find necessary and sufficient conditions for stability 
when the background magnetic field is either parallel or perpendicular to the thermal and composition gradi- 
ents that are parallel to the gravitational field. We provide approximate solutions for all the relevant modes 
involved, which are driven by gravity, conduction, and diffusion. We discuss the astrophysical implications of 
our findings for a representative galaxy cluster where helium has sedimented. 

Subject headings: galaxies: clusters: intracluster medium — instabilities — magnetohydrodynamics 
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1. INTRODUCTION 

Despite the fact that magnetic fields in galaxy clusters are 
too weak to be mechanically important, they can play a fun- 
damental role in the dynamical stability of the dilute gas by 
channelling the transport of heat. The collision-less charac- 
ter of the hot intra-cluster medium (ICM), which is generi- 
cally characterized by stable entropy gradients according to 
Schwarzschild's criterion (Piffaretti et al. 2005; Cavagnolo 
et al. 2009), enables the action of magnetic instabilities that 
are sensitive to temperature gradients (Balbus 2000, 2004). 
In particular, the magneto-thermal instability (MTI) oper- 
ates when magnetic field lines are orthogonal to a tempera- 
ture gradient parallel to the gravitational field (Balbus 2001), 
whereas the heat-flux-driven buoyancy instability (HBI) acts 
when magnetic field lines are parallel to a temperature gra- 
dient anti-parallel to the gravitational field (Quataert 2008). 
Numerical simulations suggest that the non-linear evolution 
of these instabilities leads to magnetic field configurations that 
tend to suppress the ability of the plasma to tap into the free 
energy supplied by the background temperature gradient (Par- 
rish et al. 2009; Bogdanovic et al. 2009; Ruszkowski & Oh 
2010). These magnetic field configurations can nevertheless 
support overstable gravity modes (Balbus & Reynolds 2010). 

While the landscape of thermal instabilities that render ho- 
mogeneous, dilute plasmas unstable has been well explored 
(Kunz 20 11), and even extended to account for the effects of 
cosmic rays (Chandran & Dennis 2006), very little is known 
about the effects that composition gradients can have on the 
stability of the dilute ICM. If magnetic fields do not prevent 
the efficient diffusion of ions (Narayan & Medvedev 2001; 
Chuzhoy & Nusser 2003; Chuzhoy & Loeb 2004) then the 
gradients in mean molecular weight can be as important as the 
gradients in temperature (see Section 4 and Qin & Wu 2000; 
Peng & Nagai 2009; Shtykovskiy & Gilfanov 2010; Bulbul 
et al. 201 1) and provide another source of free energy to feed 
instabilities. In order to obtain a more complete picture of the 
stability properties of the ICM, it is thus important to relax the 
assumption of a homogeneous plasma. 

As a first step towards understanding the role of composi- 
tion gradients in dilute plasmas, where magnetic fields play a 
key role by channelling the conduction of heat and the diffu- 
sion of ions, we investigate the stability of modes for which 
heat conduction is slow compared to the dynamical timescale 
involved. These modes contain, among others, the overstable 
g-modes studied by Balbus & Reynolds (2010). 



2. MODEL FOR THE MULTICOMPONENT, DILUTE ATMOSPHERE 

2. 1 . General Considerations for the Plasma Model 

In order to highlight the physical phenomena that emerge 
when composition gradients are accounted for, we focus our 
attention on a dilute binary mixture (e.g., hydrogen and he- 
lium)' in a fixed gravitational field described by- 

|^ + V-(pt.) = 0, (1) 

*_iv.(p + f,-|!ii)+., (2) 

at p \ Sir Att J 

— = Vx(^;xB), (3) 

ot 

:^|(lnP.-') = (p.-„)|l„^-V^Q.,(4) 



Here p is the mass density, v is the fluid velocity, B is the 
magnetic field, g is the gravitational acceleration, and 7 is the 
adiabatic index. The Lagrangian and Eulerian derivatives are 
related via d/dt = d/dt + v-W. 

Equations (l)-(5) describe the dynamics of a binary mix- 
ture in the low-collisionality regime and they differ from stan- 
dard MHD in three important respects: 

(/) The pressure tensor P = p±\ + — p±)bb, is 
anisotropic; where the symbols ± and || refer to the direc- 
tions perpendicular and parallel to the magnetic field, whose 

direction is given by the versor b = B / B. 

(ii) Heat flows mainly along magnetic field lines, because 
the electron mean free path is large compared to its Larmor 
radius. This process is modeled by the second term on the 

right hand side of Equation (4) via = — V)r, where 
T is the plasma temperature, assumed to be the same for ions 
and electrons, and x ~ 6 x IG^^T^/^ ergs cm^'^ s^'^ K^'^ is 
the thermal conductivity (Spitzer 1962; Braginskii 1965). 



' It is straightforward to generalize the present analysis to A'^ species. 
^ For the sake of simplicity, we do not consider relatively weaker effects, 
such as thermo-difussion, baro-difussion, etc. (Landau & Lifshitz 1959). 
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(Hi) The composition of fluid elements can change due to 

particle fluxes. Considering = —Db{b-'V)c on the right 
hand side of Equation (5) ensures that the diffusion of ions is 
mainly along magnetic field lines. This is a good approxima- 
tion when the plasma is dilute enough for the ion mean free 
path to be large compared to the ion Larmor radius. Note that 
the concentration c is related to the mean molecular weight via 
l/ji = (1 - c)(l + Zi)/^ti + c(l + Z2)/H2, where and Z^, 
with i = 1, 2, are the molecular weights and the atomic num- 
bers for the two ion species. The isotropic part of the pressure 
tensor is thus P = 2pj_/3 + p\\/3 = pk-gT / fimti, where fee 
is the Boltzmann constant and toh is the atomic mass unit. 



2.2. Background State and Equations for the Perturbations 

We assume a plane-parallel atmosphere in a constant gravi- 
tational field g = —gz which is stratified in both temperature 
and composition along the vertical direction. We consider a 
background magnetic field which is weak enough that the me- 
chanical equilibrium of the atmosphere, with scaleheight H, 
is maintained via dP/dz = —gp- In general, the background 

heat and particle fluxes do not vanish, i.e., b-'VT ^ and 

6- Vc 7^ 0, unless the magnetic field and the background gra- 
dients are orthogonal. The existence of a well defined steady 
state, i.e., V-Qg = V-Q^. = 0, demands that the background 
gradients should be linear functions of the distance along the 
direction of the magnetic field. However, even if this condi- 
tion is not strictly satisfied, the dynamics of the modes that 
we consider is unlikely to be significantly affected if the lo- 
cal dynamical timescale is short compared to the timescale in 
which the entire system evolves (see also Quataert 2008). 

The modes of interest have associated timescales that are 
long compared to the sound crossing time and it thus suffices 
to work in the Boussinesq approximation. In this limit, the 

equations for the linear perturbations b ^ ^ot+ik-x i^g^ome 



aSv = —g—z ~ ikv^^^ ' ' " 



P 



P 13 B 



SB 



2v 

a5B ~ik\\B5v , 

dp , 

a — = di'z + 

P 9 



(6) 
(7) 



7 II T 



-z- Kk- I — : — Sb;,b + bz — ; ^ I , (8) 
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dz 



5v._-Dk^^-^ 



dz B 



-iDk- ( —^5b,b + hz 



dz 



dz B 



(9) 



In agreement with the Boussinesq approximation, the velocity 
perturbations satisfy k-6v = and the fluctuations in density, 
temperature, and mean molecular weight are related via 



5p ST Sp 
p T p 



0. 



(10) 



Here, we have introduced the Alfven speed, t>A = 
the thermal speed, Wth = \/2P/p, the plasma /? 
the viscosity u of the binary mixture, the thermal diffusion 



/v 



coefficient k = xT/P, and the Brunt— Vaisala frequency 



7V2^ y^lnPp-7 = 9^1nf^^_LlV (11) 



gd_ 

7 dz 



p 



For completeness, we define here several quantities that 
play an important role in the stability analysis. We denote by 
cJdyn = and CJA = k-VA the dynamical and Alfven 

frequencies, respectively. The inverse time-scales t^^ = 
K{k-bf{-i - l)/7, = D{k-bf, and r"! ee 3fc2ut\/2i^ 

characterize, respectively, the diffusion of heat, particles, and 
momentum along magnetic field lines. We also define two 
quantities which appear naturally when thermal and composi- 
tion gradients are considered 



Nr 



£ln(r^), N^,^g:i-\n\^ 



dz 



(12) 



3. LINEAR MODE ANALYSIS 



We are concerned with modes for which heat conduction is 
slow compared to the dynamical timescale and thus cjdyn ^ 
. Because the timescales associated with viscous pro- 
cesses, which are comparable to those involved in diffusion, 
are longer than the conduction timescales by an order of mag- 



nitude (Kunz 201 1), we also have r 



> T- 



We fur- 



ther focus our attention on modes for which magnetic tension 
is unimportant, i.e., fc||_ff <C and thus the Alfven fre- 

quency is small compared to other inverse timescales. Since 
the plasma (3 ranges from lO'^ from the centers of cool core 
clusters to 10* in the outskirts of the ICM (Carilli & Taylor 
2002), there is a reasonable range of wavenumbers for which 
a local analysis is sensible. We can thus summarize the regime 
in which we are interested according to 



Wdyn > 



(13) 



For a homogeneous plasma, the modes satisfying these condi- 
tions encompass the overstable g-modes studied in Balbus & 
Reynolds (2010). The set of Equations (7)-(10) allows us to 
address the behavior of these, as well as other new modes, in 
the presence of a non-vanishing gradient in the mean molec- 
ular weight and account self-consistently for the diffusion of 
ions along magnetic field lines. 

3.1. Heat Diffusion Along Vertical Magnetic Field Lines with 
No-Ion Diffusion 

Let us first consider the case in which there is no ion dif- 
fusion. Setting D = I' = the dispersion relation corre- 
sponding to Equations (7) to (9) factorizes and possesses as 
non-trivial solutions two Alfven waves, with a = ± jwa, and 
the roots of the polynomial 



+ aia^ + a2<y + 03 = , 



with coefficients 

-1 



ai 



k^ ' 



0.3 



(14) 

(15) 
(16) 



where the approximate expressions for the coefficients follow 
from the general considerations outlined above. 
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(a) gravity modes (b) conduction modes (c) diffusion modes (d) gravity modes (e) conduction modes 




B = B^:D=0 B = B,:D = Q B = B,:D^O B = B^ B = B^ 



Fig. 1. — Graphic representation of the stability of modes for which conduction is slow compared to the dynamical timescale, i.e., uj^yn 3> . The 
various panels show the unstable regions (gray) for each of the modes that can be excited when the background magnetic field is parallel (a, b, and c) or 
perpendicular (d and e) to the background thermal and composition gradients. The solid line corresponds to A'^'^ = 0; the horizontal dotted line represents 
d IriT/d In P = {"f — l)/2j; and the dashed lines correspond to din T/d\n P = ±d \n^/d\nP. For B = Bz, gravity modes (a) can become either stable or 
overstable, while modes driven by conduction can become unstable (b). If ions can diffuse efficiently along magnetic field lines, a new type of mode can become 
unstable (c). For B = Bx, both gravity modes (d) and conduction modes can become unstable (e); while ion diffusion only leads to decaying modes. 



The Routh-Hurwitz stability criteria that predict exclu- 
sively negative real parts for the roots of the cubic Equa- 
tion (14) are: ai > 0,03 > 0, and aia2 — 03 > 0. The 
first condition is trivially satisfied while the other two imply 



iV|^<0, 



(18) 
(19) 



respectively. To leading order, two of the solutions of Equa- 

1/2 

tion (14) are given by ct « ±102 + (03 — aia2)/2a2, i.e.. 



k 2tc 



iV2 



(20) 



which correspond to gravity modes. In the absence of a gra- 
dient in the mean molecular weight, these reduce to the 17- 
modes discussed in Balbus & Reynolds (2010). The third root 
is given by cr « —0,3/0,2, i.e.. 



7V2 



(21) 



and corresponds to a mode driven by conduction. Assuming 
that A^^ > 0, g-modes are overstable if the condition (19) is 
not satisfied, while conduction modes are unstable if (18) is 
not fulfilled. 

In order to understand how these modes behave in the 
parameter space spanned by temperature and composition 
gradients, it is convenient to use that dlnP/dz = —1/H 
and work with the dimensionless variables d\nT/dln P and 
dhifi/dhiP. The classical requirement for stability against 
buoyancy, i.e., N'^ > becomes 



dlnT dlnfi 
d\nP d\nP 



7 - 



1 



7 



while the conditions (18) and (19) become, respectively. 



dlnT 

d\nP 
dlnT 

d\nP 



> 



< 



dln^ 
~dlnP 
7-1 

27 ■ 



(22) 

(23) 
(24) 



The first panel (a) in Figure 1 shows the regions of parameter 
space where g-modes in Equation (20) are stable, over-stable, 
or unstable (gray area). The second panel (b) shows that the 
modes in Equation (21), which are driven by conduction, can 



be either stable or unstable. Note that in the region of param- 
eter space where both gravity and conduction modes are over- 
stable/unstable they both grow with comparable rates. Panel 
(a) in Figure 2 shows the stable region (white) of parameter 
space satisfying the conditions (23) and (24) simultaneously. 



3.2. Heat and Ion Diffusion Along Vertical Magnetic Field 

Lines 

Let us now consider the situation in which ions diffuse 
mainly along magnetic field lines. If _D 7^ 0, the dispersion 
relation corresponding to Equations (7) to (9) yields as non- 
trivial roots a = ±iujA and the solutions to 



with coefficients 



= 0. 



(25) 



bi 
b2- 

bi 



',, ^ 'c ' 



fc2 ^ 



fci 

fc2 



/r2 

k^ 



(26) 



/r2 

rfc2 



Tfj, - 



(27) 



fc2 



The approximations on each of the last terms on the right hand 
side hold under the same considerations that lead to the sim- 
plified coefficients in Section 3.1. The approximate expres- 
sions for the coefficients bi and a; are identical for i — 1,2,3. 

The Routh-Hurwitz stability criteria for a quartic polyno- 
mial with real coefficients requires: 61 > 0, 64 > 0, 6162 — 
63 > 0, and - bjbi - 6| > 0. While 61 > is trivially 

satisfied, 64 > implies: 



(29) 



The other two conditions lead again to inequalities (18) and 
(19). Therefore, the effects of diffusion require that only the 
additional condition (29) be met for the system to be stable. 

To leading order, three of the roots of the dispersion relation 
(25) are given by Equations (20) and (21). The fourth solution 
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Fig. 2. — Panels (a), (6). and (c) show a graphic representation of the stable regions (white) satisfying all the Routh-Hurwitz stability criteria discussed in 
Section 3. Panel (d) is a schematic representation of the correspondence between the regions of a representative galaxy cluster with radial temperature and mean 
molecular weight profiles as shown in Figure 3 and the plane spanned by {dlnT / din P, din fi/ din P) . 



consists of a mode driven by ion diffusion 



(30) 



This is an unstable mode if either of the criteria (18) or (29) is 
unfulfilled. 

Thus, if ions can diffuse along magnetic field lines, in ad- 
dition to requiring that the gradients in temperature, pressure, 
and mean molecular weight satisfy the inequalities (23) and 
(24), we must have 



dlnT din fi 
d\nP^ d\nP 



(31) 



The third panel (c) in Figure 1 shows that the modes in 
Equation (30), which are driven by diffusion, can be either 
stable or unstable (gray). Their growth rates are estimated 
to be an order of magnitude smaller than either g-modes or 
conduction modes. The importance of these diffusion modes 
resides in that they can become unstable in regions of param- 
eter space which are stable against g-modes and conduction 
modes. Panel (b) in Figure 2 shows the stable region (white) 
of parameter space satisfying all the conditions (23), (24), and 
(31) simultaneously. 

3.3. Heat and Ion Diffusion Along Horizontal Magnetic Field 

Lines 

If the background magnetic field is perpendicular to the 
thermal and composition gradients, the dispersion relation 
that governs the stability of the atmosphere is 

cr^ + Cicr** + C2a^ + Cgcr^ + €4(7 + C5 , (32) 

where the coefficients 



fr2 

— 1 —1 I 

Cl ~ T-c + 



-1-1^ 
'c 'v ^2 ' 



(33) 
(34) 
(35) 
(36) 



' In the absence of diffusion, the only result that is modified in this section 
is that the root a = —t7 ^ for D ^ becomes a = for D = 0. 



(37) 



are subject to the same considerations employed in deriving 
the approximate expressions for the coefficients Oj and hi. 

The Routh-Hurwitz stability criteria for a fifth degree poly- 
nomial requires ci > 0, C1C2 — C3 > 0, C1C2C3 



C1C5 > 0, C4(ciC2C3 



C?C4 



C?C4 



C3) > C5(ciC^ - 2ciC4 



C2C3 + C5), and C5 > 0. In the limit under consideration, i.e., 
Wdyn ^ > Tv^^ these inequalities become, respectively. 



T/t, 



Nr 



2 



{N 
2 



>0, 



Nr 



T/fJ. 



> 0. 



(38) 
(39) 
(40) 
(41) 
(42) 



The first of these conditions is trivially satisfied. This is also 
the case for inequality (39), since it can be written as [(7 — 
I) / -f][dP / dz]'^ / P p > 0. Therefore, the only independent 
condition required for stability is -/V|,y^ > 0. 

Two of the approximate solutions to the dispersion relation 

1/2 

(32) are given by a w ±ic2 + (C3 — ciC2)/2c2, i.e.. 



fc2 I J.2 



1 

2i> 



7V2 



(43) 



which correspond to stable gravity modes. The other two so- 
lutions correspond to a conduction and a viscous (decaying) 
mode, which are, respectively 



Ar2 



1.2 



(44) 



We conclude that when the magnetic field is perpendicular 
to the temperature and composition gradients, the stability of 
g-modes requires only that A^^ > g, whereas the stability of 
conduction modes requires also 



d In T d In /i 
dlnP ^ dlnP 



(45) 



The fourth panel (d) in Figure 1 shows the regions of pa- 
rameter space where g-modes in Equation (43) are stable or 
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unstable (gray). The fifth panel (e) shows that the modes in 
Equation (44), which are driven by conduction, can be either 
stable or unstable. These regions are significantly different 
from the corresponding regions in panel (b), for which the 
direction of the background magnetic field is parallel to the 
direction of the temperature and composition gradients. Note 
that, unlike the case where B = B^, 5-modes and conduc- 
tion modes cannot be simultaneously unstable when B = B^- 
Panel (c) in Figure 2 shows the stable region (white) of pa- 
rameter space satisfying the condition (45). 

4. ASTROPHYSICAL IMPLICATIONS 

We have provided a detailed description of the various in- 
stabilities that can be present in the parameter space spanned 
by (d In T/d In P, din In P) without imposing restric- 
tions on the relative signs of the gradients involved. We can 
now frame our results, summarized in Figure 2, in the con- 
text provided by observations and theoretical models address- 
ing the temperature and composition structure of galaxy clus- 
ters. Figure 3 shows a schematic representation of the tem- 
perature and mean molecular weight profiles of a representa- 
tive galaxy cluster. The temperature profile sketched there re- 
sembles the results obtained by observations (Vikhlinin et al. 
2006), whereas the mean molecular weight profile is akin to 
helium sedimentation models (Bulbul et al. 201 1). If the peak 
in the temperature takes place at a larger radius than the peak 
in the mean molecular weight, then there are three distinct re- 
gions defined by the signs of the temperature and composition 
gradients. These regions correspond to different quadrants in 
the (din r/dlnP, din /i/dln P) plane as shown in panel (d) 
in Figure 2. 

The joint analysis of Figures 1 and 2 allows us to under- 
stand the implications that mean molecular weight gradients 
can have for the various regions of a representative galaxy 
cluster as depicted in Figure 3. If the magnetic field in the in- 
ner ICM is perpendicular to the background temperature gra- 
dient dT/dz > 0, as suggested by the end states of initial 
configurations with B = Bz that are HBI-unstable but evolve 
to B ^ Bx, then this region is stable if d/i/dz = (although 
there could be overstable g-modes driven by radiative cooling 
Balbus & Reynolds 2010). However, if df^i/dz > this region 
is unstable to conduction modes. If the magnetic field in the 
outer ICM is parallel to the temperature gradient dT/ dz < 0, 
as suggested by the end states of initial configurations with 
B = Bx that are MTI-unstable but evolve to B w Bz, then 
this region is overstable to g-modes and unstable to conduc- 
tion modes even if d/x/dz = 0. A negative gradient in the 



mean molecular weight is unable to stabilize these modes and 
it can further drive unstable diffusion modes. Therefore, even 
if the outer/inner ICM is stable to both the MTI and the HBI, 
it could be rendered unstable to modes driven by either grav- 
ity, conduction, or diffusion. Unless the magnetic field is or- 
thogonal to the background gradients, the region that we have 
denoted as intermediate is likely to be unstable to conduc- 
tion and diffusion modes, but not to gravity modes. These 
g-modes, however, could be relevant if the temperature pro- 
file peaks at a smaller radius than the mean molecular weight 
profile. 

We will address the limit in which conduction is fast com- 
pared to the dynamical timescale in a forthcoming paper. This 
limit contains the generalization of the instabilities that be- 
come the MTI and the HBI in the limit of a homogeneous 
plasma. 
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Fig. 3. — Schematic representation of the mean molecular weight (red solid 
line) and temperature (black dashed line) profiles of a representative galaxy 
cluster as suggested by observations (Vikhlinin et al. 2006) and theoretical 
models (Bulbul et al. 2011). The regions denoted by "inner", "intermediate", 
and "outer" ICM (delimited by dotted gray lines) correspond to three differ- 
ent quarters in the (d In T/dln P, d In/i/d In P) plane (Fig. 2.). The mean 
molecular weight for a homogeneous cluster with primordial abundance is 
fi ^ 0.59. 
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